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ABSTRACT 


The implementation of Satellite based Doppler positioning systems fre- 
quently requires the recovery of transmitter position from a single pass of 
Doppler data. The least squares approach to the problem yields conjugate solu- 
tions on either side of the satellite subtrack. It is important to develop a pro- 
cedure for choosing the proper solution which is correct in a high percentage 
of cases. 

This paper derives a test for ambiguity resolution which is the most power- 
ful in the sense that it maximizes the probability of a correct decision. When 
systematic error sources are properly included in the least squares reduction 
process to yield an optimal solution, the test reduces to choosing the solution 
which provides the smaller valuation of the least squares loss function. When 
systematic error sources are ignored in the least squares reduction, the most 
powerful test is a quadratic form comparison with the weighting matrix of the 
quadratic form obtained by computing the pseudo-inverse of a reduced rank 
square matrix. 

The paper provides a formula for computing the power of the most power- 
ful test. A numerical example is included in which the power of the test is com- 
puted for a situation which may occur during an actual satellite aided search 
and rescue mission. 
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INTRODUCTION 


A Search and Rescue Orbiting Satellite System has been suggested as a 
means for locating distress signals from downed aircraft, small boats, over- 
land expeditions, etc. 1,2,3> Positioning with this system would be accomplished 
using emissions from Emergency Locator Transmitters (ELT), now available on 
most U.S. aircraft. An essential requirement of the system is that the location 
of ELT emissions be determined to within a few kilometers using data obtained 
during a single pass of a 700-1200 kilometer altitude satellite. 

One specific mission proposal 5 involves the incorporation of search and 
rescue mission instruments on a TIROS-N satellite in a circular, 850 Km alti- 
tude, polar orbit. The instruments would be capable of detecting and locating 
ELT’s operating at 121.5 and 243 MHz, as well as new and improved ELT's 
operating on the 406 MHz frequency authorized for ground to satellite use by the 
1971 World Administrative Radio Conference on Space Telecommunications, me 
nominal launch data is J;in. 1981. 

The requirements of such a mission lead to a certain estimation problem: 
that of determining transmitter position from a single pass of Doppler data and 
a satellite epoch state within or near the data span. In fact, it is very difficult 
to recover three position components from a single pass of Doppler observations. 
The usual procedure is to assume the transmitter to be on the surface of an 
ellipsoidal Earth and to use a conventional least squares filter to estimate longi- 
tude and latitude components of position. Also, the carrier frequency cannot be 
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assumed as known. Hence it is necessary to recover the carrier frequency from 
the data along with transmitter longitude and latitude. 

All least squares reduction processes are equivalent to finding a local min- 
inum of a surface. In the case of single pass Doppler position determination the 
multidimensional surface in question has two local minima, representing solu- 
tions on each side of the satellite subtrack. Hence two solutions are available 
and they can be obtained by choosing initial values for the least squares itera- 
tion process on either side of the satellite subtrack. Within the context of a 
satellite aided search and rescue mission the ability to choose the correct least 
square solution can be critical. Unless one develops a decision procedure which 
is correct in a high percentage of cases, situations may occur in which it will be 
necessary to send rescue parties to two widely separated locations thus increas- 
ing the cost of each rescue. The ambiguity problem, the distinguishing of the cor- 
rect from the spurious solution, is the subject of this paper. 

The problem arises in other applications. An example is the proposed satel- 
lite-aided coastal surveillance in support of the recently imposed 200 mile fishing 
limit. 6 This system would require a single pass Doppler positioning capability. 
Also, Green 7 has analyzed the ambiguity problem in the course of a comprehen- 
sive treatment of the NIMBUS-TWERLE experiment. But he does not attempt 
to derive the optimal decision procedure. In addition, the technique recommended 
for computing the probability of a correct decision which relies on a cumulative F 
distribution is incorrect. 
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In what follows it is shown that the residual sets of the two solutions have 


different statistical properties and that this fact can be exploited to develop a 
simple quadratic form test which is the most powerful in the sense that it max- 
imizes the probability of a correct decision. 

A numerical example is given in which the power of the test is computed for 
a situation which may occur during an actual satellite aided search and rescue 
mission. 

RESIDUAL STATISTICS 

In this section we develop the statistical properties of the a-posteriori least 
squares residuals for both the correct and the spurious solution of the single pass 
Doppler positioning problem. The treatment is sufficiently general to account for 
systematic error sources which are left unadjusted in the reduction process. 

Let x be the three dimensional vector representing true values of transmitter 
longitude, latitude and transmitter carrier frequency. Let y be an m dimensional 
set of correct or noiseless Doppler obseiwations obtained from a single satellite 
pass. The actual observations are corrupted by noise. Hence 

y = y + v, K(v) = 0. E(w T ) = Q (1) 

where Q is a full rank covariance matrix, and v is normally distributed. In addi- 
tion, there are other parameters z (orbit epoch state, carrier frequency drift rate, 
ete. s ) which influence the functional relationship between y and x. If tiie estimates 
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of these parameters are uncertain to a significant extent, this fact must be taken 
into account in obtaining an optimal estimate x of x. Let z be an estimate of 2. 
Then 

z = z + a. E(a) = 0, E(aa T ) = P (2) 

where P is full rank and i normally distributed. The functional relationship be- 
tween y, X and 2 can be expressed as 

y = f (x, z ) (3) 

Let x n be a nominal estimate of x. Expand equation 3 in a first order Taylor 
series about x^, and z and include the results of equations 2 and 1. The result is 

Sy = ASx + y (4) 

where 


7-Ba+y,E(y) = 0, E(y> T ) - BPB T + Q = 0 

(5a) 

8y = y - f . z) 

(5b) 

bx = x - x., 
N 

(5c) 

A=— , 

B.” 

(5d) 

Hi 
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Equation 4 can be used as an equation of condition for a standard least squares 
estimator. The result is 

x = x N + ( A T ^" , A) -1 A T 0‘ 1 Sy (6) 

The left side of equation <3 can be used as a new nominal and the process can be 
repeated until a convergence criterion is satisfied. By employing the usual 
linearity assumptions one can show that the left side of eq. 6 is an optimal esti- 
mator of x in the sense that the trace of its covariance is a minimum. 9 
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We require that the first and second order statistics for the a-posteriori 
residuals of the least squares filter when the iteration process converges to a 
correct solution. Assume that the nominal value x used as an initial guess 
for the least squares iteration process is sufficiently close to £ that the first 
order Taylor series expansion represented by eq. 4 is valid. Then by combining 
eq. 4 with eq. 6 we obtain 

x = x + (A t <£ _1 A ) -1 A T 0'* y (7) 

The a posteriori residuals are defined as 

R = y - f ( x, z ) ( 8 ) 

Assume that x is a good estimate of x. and rewrite equation 8 as 

R = f (x, z ) - f (x, z) + 7 ( 9 ) 

Then a first order Taylor series expansion about x and 2 yields 

R = [I - A(A T ;- 1 A)" I A T ^- 1 ]7 (10) 

Hence It is a linear function of a normal random variable and, 

E(R)= 0 (Ha) 

COV(R) = E(RR T ) = <p - ACA^^A^A^ C (lib) 

In fact the covariance matrix C of R is of reduced rank. This can be seen by 
focusing attention on eq. 10. Because of the interpretation of the least squares 
reduction process as a projection operator 11 it follows that the matrix which 
maps > onto R is a projection operator on the eomplement space of the range 
of A. Since y is a full rank random vector and since the range space of A is of 
rank 3, it follows that the M dimensional covariance matrix C is of reduced rank 
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IV1-3. These facts can be summarized in the following 


Theorem: 

Assume that a least squares process has been implemented to solve the 
single pass Doppler positioning problem and that the process has converged to 
a point near the true transmitter position. Assume also that the equation of 
condition (eq. 4) properly accounts for systematic error sources which are left 
unadjusted in the least squares process. Then equation 6 provides a minimum 
variance estimate of transmitter posit. on. The a posteriori residuals of the 
process have first and second order moments given by eq. 11a. , and lib. The 
covariance matrix of the residuals is of reduced rank M-3 where M is the num- 
ber of data points. 


To obtain residual statistics for the spurious solution, let x' be the expected 
value of the least squares solution when the iterative process converges to the 
wrong side of the satellite subtrack. The spurious solution is obtained when the 
initial estimate x N is within a certain radius of convergence of x'. Hence eq. 6 
can be considered as the least squares estimator for the spurious solution but 
with the partial derivative matrices evaluated at x' rather than x . Perform the 
following first order Taylor series expansion 

f(x N .z) - ffx'.'zj = A(x N - x' ) + B(z , S') (12) 

where the partial derivative matrices A and B are evaluated at x' . When eq. 12 
is combined with eq. G, the result is 
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x = x' + ( A t »£ -1 \)- 1 A t 4>~ 1 [ f (x,z) - f ( x' , 1? ) + ' ^ (13) 

But since x' is interpreted as the expected value of the least squares estimator, 
we have 

x = x' + (A T vp“ 1 A)" 1 A (14) 
By combining eq. 9 with eq. 14 and employing another Taylor series expansion 
about x ' , we obtain an expression for the residuals of the spurious solution as 
R s = f(x,z ) - f(x\z ) + [I - A( A t <£ -1 A ) -1 A t sP _1 ] y (15) 

Hence the first and second order statistics for the a-posteriori residuals of the 
spurious solution are 


E(R S ) = f (x ,z) - f ( x',7 ) = U 

(16a) 

COV(R s ) = <p - A( A T 0 -1 A) -1 A T 

(16b) 


The symbol U is the difference between the true Doppler signal received from a 
transmitter located at its actual position represented by x and the two Doppler 
signals which would be received if the transmitter were at the position repre- 
sented by x' . 

The covariance matrix of the residuals for the spurious solution given by 
16b has the same form as the covariance matrix for the residuals of the true 
solution as given by lib. The difference is that the partial derivative matrices 
A and B defined in equation 5b. are evaluated at approximately mirror image 
points relative to the satellite subtrack. However, numerous simulations have 
shown that the matrix form given by equation lib or 16b are nearly identical 
when evaluated at mirror image points. This is the result of the near symmetry 
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of the function defined by equation 3 about the satellite subtrack. These facts 
are summarized as follows: 


Theorem: 

The covariance matrix of the a-posterlori residuals of the spurious solul.on 
is the same as that of the a-posteriori residuals for the correct solution. The 
expected value of the residuals associated with the spurious solution is the dif- 
ference between the signal received if the transmitter were at x and the signal 
received if the transmitter were at x' . 


MOST POWERFUL TEST TOR AMBIGUITY RESOLUTION 
Assume that a least squares solution to the single pass Doppler positioning 
problem has been implemented and that solutions on sides a and b of the satellite 
subtrack have been obtained. Let r and r. be the residual sets of the two solu- 

H 0 

tions. We require an algorithm which permits us to decide which of the two 
residual sets is associated with the correct solution. Furthermore, it is to be 
the most powerful test in the sense that among all possible such tests, it maxi- 
mizes the probability of a correct decision. 

The problem can be formulated in terms of a standard testing of a statistical 
hypothesis. Let the null hypothesis be that the vector r t is a realization of a 
random vector R whose mean is 0 and whose covariance matrix is C as defined 
by eq. lib. The null hypothesis is equivalent to the assumption that the correct 
solution is on side a. The alternative hypothesis is that r n is a realization of 
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a random vector R s whose expected value Is U as defined by equation 16a. and 
whose covariance matrix is C. The resultant test is to be the most powerful in 
the sense that it maximizes the probability of rejectingthe null hypothesis when 
in fact the null hypothesis is false. 

The derivation of the most powerful test would be quite straightforward if 
it were not for the fact that the random vectors R and 1^ are of reduced rank. 

It is necessary to account for thi3 fact by performing a linear transformation 
of variables by means of a matrix T of dimension (M-3) by M. The transforma- 
tion defined by T must have the property that the composite mappirg TC have a 
range space of dimension M-3. (One such transformation is defined by a matrix 
whose rows are the M-3 eigenvectors associated with the non-zero eigenvalues 
of C.) Define 


£ = TR 

(17a) 

£ S =TR S 

(17b) 

V « = T V V b = Tr b 

(17c) 

W 2 TU 

(17d) 


The null hypothesis is now that the vector V o is a sample of the M-3 dimen- 
sional random vector £ ’/hose mean is 0 and whose covariance is 

COV(£ ) = TCT t = D (18) 

where D is a full rank matrix. The alternative hypothesis is that V b is a sample 
from the random vector £_ with mean W and covariance ! >. 


Imagine an experiment whose first step is to choose either random variable 
r or random variable . s with probability .5 for each choice. Next obtain a real- 
ization V from whichever random vector was chosen at the first 3tep. The prob- 
ability structure of the experiment can be described as follows: Define a two 
point random variable x as x = 0 if random vector.' is chosen and x = 1 if ran- 
dom vector . s is chosen. The outcome of the expe. *ment can be described as a 
two tuple (V, x) where the value of x describes the chosen probabil.iy law and 
where V is the realization. From equations 10, 15 and 17a, 17b, it is seen that 
random vectors and £ s are linear functions of normal random vectors. Since 
they are also full rank they are normally distributed. Hence the joint probability 
density function of the experiment outcome is 


i ( V , x ) = 


1 exp* ^ Iv t d ' 1 v] 




(M-31,'2 D 1 2 


(19) 


if x - 0 


1 ( V . x ) = 


J CX p- ■{ UV-W^T )* 1 (V -*)] 


2 d 


1 2 


if x = 1. 


The Neyman- Pearson Lemma 10 asserts that in the case of simple alternative 

hypotheses, a maximum likelihood criterion is the most power- 
ful test for accepting or rejecting the null hypothesis. An 
application of Bayes' theorem** yields the conditional probability 
density function of x given a realization V a of V as 
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f ( x I V = V a ) = 


exp" 2 [vJd-W,] 

exp" 2 [v>" lv J+exp _ 2 [ < v a - w) T D _1 (v a - w>] 


( 20 ) 


if x = 0 


f(x I V = V a ) = 


exp 


f 2 [(V a -W) T D _1 (V a -W)] 


-I [v^D -1 V a ] + exp - £ [(V a - W) T D* 1 (V a -W)] 


exp 


if x = 1. 


From equation 20 it follows that the maximum likelihood criterion is 


f (Ol V = V a ) > f ( 1 1 V = V a ) 


(21) 


But since we have alternative and mutually exclusive hypothesis: 


f ( 1 I V = V a ) = f (0 I V = V b ) 


( 22 ) 


From equations 20, 21, and 22 we obtain the result that the maximum likelihood 
test is to accept the null hypothesis if and only if 

VjD- 1 V a <VJD-i V b (23) 

Equation 23 can be written as 

r T t t d -i t ^ < r t t t d-*T r b (24) 

But in fact, 

T T D- 1 T = C t (25) 

where C denotes the unique pseudoinverse 9 of C, This can be seen by showing 
that the left side of equation 25 satisfies the four axioms of the pseudoinverse. 

cc f c = c (26a) 

C’ CC f = C’ (26b) 
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(ccV = ex’ 


(C T C) T = C T C 


Hence the criterion for accepting the null hypothesis is 


r T C r r < r. T C T r. 

a o b b 


A further simplification is possible. 


Theorem: 


Both residual set r a and residual set r b are in the range space of C. 


Proof: 

Residua 1 sets r a and r b are realizations of random vectors R and R s which 
are defined by equations 10 and 15. Hence it is sufficient to show the existence 
of random . uriables G and G s such that 

R = CG (28a) 

R S = CG S (28b) 

From equations 10 and lib, it is clear that the required representation for G is 
0' 1 7. Also notice that from equations 13 and 14 

(AV 1 A) -1 A t 0 _1 [f(x,z) - f( x', z' )J = 0 (29) 

when x’ is the expected value of the least squares estimator when the iteration 
process converges to the wrong side of the satellite subtrack. Hence equation 15 
can be rewritten 


R s = [I - A( A 1 *" 1 A ) -1 A t £ _1 ] [ f ( x,7 ) - f ( x\z ) + y] 
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The correct expression for G s is 

G s = <£[f (x,?) - f( k',z ) + y] (31) 

Theorem: 

Let C be a matrix transformation whose form is that of lib, and let X be a 
vector in the range space of C. Then 

X T C t X=X T 0- 1 X (32) 

Proof: 

Let F be such that 

X = CF (33) 

where C has the form of equation lib. The following set of equations proves the 
result 


X T CX = F T CC t CF 

(34a) 

F T cc t CF = F T CF 

(34b) 

F t CF = F J Ccp- l CF 

(34c) 

F T C<p~ i CF = X T 0" l X 

(34d) 


The above results imply that inequality 27 is satisfied if and only if 

r J • > ' 1 r a < r b T< ^‘ 1 r b ( 35 ) 

We have derived the fundamental result of the paper, 

Theorem: 

Assume that the least squares solution to the single pass Doppler positioning 
problem has been implemented and that solutions on either side of the satellite 
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subtrack have been obtained. Let r a and r () be the residual sets of the two solu- 
tions. Assume also that systematic error sources have been properly modeled 
in the reduction process so that when ambiguity is properly resolved the result- 
ing least squares estimator is optimal in a minimum variance sense. Then the 
most powerful test for ambiguity resolution is obtained by choosing the solution 
on side a of the satellite subtrack as the correct transmitter position if and only 
if inequality 35 is satisfied. The quadratic forms on the left and right sides are 
the valuations of the least squares loss function obtained at each of the two solu- 
tions and, hence, are available as an automatic byproduct of the least squares 
reduction process. If no systematic error sources are present and if the data 
noise is uncorrelated and of constant variance, the most powerful test reduces 
to choosing the solution whose residuals provide the smaller square sum. Once 
the residuals are available the probability of the test providing a correct answer 
can be obtained by computing 

exp' i V 

P = (36) 

exp" 5 ' R J + exp ’ 3 R i>° R t. 

where the required probability is P if P > .5 and 1-P if P < .5. 

An entirely different approach which also derives the above result is given 
in the appendix. 

If the systematic error sources defined by equation 2 are neglected in the 
least squares reduction process, the residual statistics are different from those 
given by equations 11a, lib and 16a, and 16b. Hence for this situation the above 

I 
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theorem is invalid. The most powerful test for this case can be obtained by 
noticing that if the uncertainties implied by equation 2 are ignored in the least 
squares reduction process, the least squares estimator becomes 

x = x N + (A T Q"*A) _1 A T Q -1 6y (37) 

By the same reasoning that obtained equation 10 from equation, 6, an equation for 
the residuals of the correct solution can be obtained as 

R = [I . A(A t Q- 1 A)- 1 A t Q-» )y (38) 

Define 


H = I - A(A t Q -1 A) -1 A t Q -1 (39) 

Then residual statistics for the correct solution become 


E(R) = 0 (40a) 

COV(R) = H^H t = C 0 (40b) 

By repeating the logic which obtained equation 15 from equation 6 we obtain an 
expression for the residuals of the spurious solution 

R s = f (x\ z) - f ( x',7 ) + Hy (41) 

Hence 

E(R S ) = f (x , z) - f (x',7 ) = U (42a) 

COV(R s ) = H£H t = C 0 (42b) 

The reasoning employed from equation 17 to equation 27 is again applicable and 
we have the following result, 
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Theorem: 


Assume that the least squares solution to the single pass Doppler position- 
ing problem has been implemented and that solutions on either side of the satel- 
lite subtrack have been obtained. Let r a and r b be the residual sets of the two 
solutions. Assume that the uncertainty associated with parameters defined by 
equation 2 have been ignored in the least squares reduction process. Then the 
most powerful test for ambiguity resolution is obtained by choosing the solution 
on side a of the satellite subtrack as the correct transmitter position if and only 


if 

r T Cl r < r.C« r. 

a 0 a b 0 b 

where C Q is defined by equation 39b. 


(43) 


It would be inconsistant to ignore systematic error sources in a least squares 
data reduction and then include their effects in an effort to resolve ambiguity. 
Hence it is unlikely that the test implied by inequality 43 would be implemented in 
practice. However, the test is of use in providing a bound on the probability of 
ambiguity resolution in the presence of different types and magnitudes of sys- 
tematic error sources. 

COMPUTING THE POWER OF THE TEST 
Suppose that an optimal least squares estimator defined by equation 6 has 
been implemented to determine transmitter position from a single pass of Dop- 
pler data. The inequality 35 provides the most powerful test for ambiguity 
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resolution. The power of the test is the probability of inequality 35 being satis- 
fied when the correct solution is on side a. 

From equations 10 and 15 it follows that 

iU=r o -r 6 (44) 

With regard to determining the power of the test it will be seen that the sign 
ambiguity is irrelevant. We will assume the positive sign. Inequality 35 can 
now be rewritten as 

r > _1 r fl < (r a -U)T*-i(r a -U) (45) 

Expand inequality 45 to obtain 

2U T <p~ l r a < U T <p~ l U (46) 

If the correct solution is on side a, the probability of satisfying inequality 45 is 
the probability that a standard normal random variable satisfies inequality 


„ u V 1 u 

x < 

2(U t 0" 1 C0- 1 U) 1/2 


(47) 


where C is defined by equation lib. But from equation 29, and lib it follows that 

C0 _1 U = U (48) 

Hence inequality 47 is equivalent to 


< CU^V 1 !]) 1 2 
2 


(49) 


The power of the test can be obtained by computing the right side of inequality 49 
and consulting a tabulation of the cumulative normal distribution for the proba- 
bility of the standard normal variable not exceeding the number in question. If 
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no systematic error sources are present and if the data noise is uncorrelated 
and of constant variance, inequality 49 reduces to 


< 

2 a 


1/2 


(50) 


where a is the standard deviation of the data noise. 


NUMERICAL EXAMPLE 


To illustrate the application of equation (50), an example that is typical of 
the satellite-orbit search and rescue system to be tested using the TIROS-N 
satellite 1 - 2 was chosen. A circular satellite orbit with a radius of 7200 
kilometers and an inclination of 98.7° was simulated. The latitude of the 


transmitter was about 36°N, and the distance from the satellite subtrack 
was chosen so that the satellite elevation was 20° at the time of closest approach. 
The range rates that would be measured in such a configuration using one min- 
ute intervals between observations over a span of thi’ee minutes were calculated, 
and are listed in the second column of both Table I and Table II. The first column 


in each of these tables, the time of observation, is given relative to the time of 
closest approach of the satellite to the transmitter, which occurs at the second 
range-rate observation. 

In Table I, noise was generated using a standard deviation of about 2.9 meters 
per second, which, at a transmitter frequency of 40GMHz, corresponds to a de- 
termination of the received frequency to within a bin width of 10 Hz. The noise 
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values obtained are listed in the second column of Table I. The simulated range 
rate observations were formed by adding the range rates in the first column to 
the corresponding noise values in the second. Using these simulated observations 
a least squares position location algorithm was used to solve both for the mini- 
mum near the transmitter position, and also for the false minimum on the op- 
posite side of the satellite subtrack. The residuals of the least squares fits to 
these minima are listed in the fifth and fourth columns of Table I respectively. 
The vector U appearing in equation (50), which is defined in equation (44) as the 
difference between the residuals at the false minimum and those at the minimum 
near the true position, is listed in the last column. The fundamental test, equa- 
tion (35), in this case gives the correct results since the sum of one squares of 
the fourth column, 23.684, exceeds that of the fifth column, 9.027. Using the 
residual difference in equation (50), the right side of (50) is equal to 0.438 
which, consulting a table of the cumulative normal distribution, yields a value 
of 67% for the power of the test, a low value. 

Tn Table II the process above has been repeated, except that in this case the 
receiver resolution is taken to be 1 Hz at 406 MHz. The power of the test in- 
creases to 99.99%. Clearly the power is a strong function of data noise standard 
deviation. It is also a function of Beacon latitude, satellite elevation at TCA, and 
the data distribution along the satellite pass. 
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Summary 


The use of Satellite Based, single pass Doppler positioning systems creates 
a certain estimation problem - that of determining transmitter position from a 
single pass of Doppler data and a satellite epoch state within or near the data 
span. A least squares approach is usually employed to solve the problem. How- 
ever, all least squares reduction processes are equivalent to finding a local min- 
imum of a surface. In the case of single pass Doppler position determination the 
multidimensional surface in question has two minima, representing solutions on 
each side of the satellite subtrack. In many applications it is important to have 
a procedure for choosing the proper solution which is correct in a high percent- 
age of cases. 

It is shown that the residual sets of the correct and the spurious solutions 
have different statistical properties and that this fact can be exploited to develop 
a simple quadratic form test for ambiguity resolution. This test is the most 
powerful in the sense that it maximizes the probability of a correct decision. 
When systematic error sources are properly included in the data reduction 
process to yield an optimal solution, the test reduces to choosing the solution 
which provides the smaller valuation of the least squares loss function. When 
systematic error sources are ignored in the least squares reduction, the most 
powerful test is a quadratic form comparison with the weighting matrix of the 
quadratic form obtained by computing the pseudo-inverse of a reduced rank 
square matrix. 


20 


The paper provides a formula for computing the power of the most powerful 
test. The power of the test is a function of the data noise, the number and dis- 
tribution of the data points, transmitter latitude, and the maximum elevation of 
the data points, transmitter latitude, and the maximum elevation of the satellite 
pass relative to the transmitter position. A numerical example is Included in 
which the power of the test is computed for a situation which could occur during 
* an actual satellite aided search and rescue mission. 
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Appendix 


Alternate Derivation 

The probability that the true position x of the transmitter lies in a given 
region A when a particular series y of range rates are observed is given by 
[13] 

P[x in A l y ] = — — - f p( y| x ) p( x )dx (Al) 

p( y> J A 

where p(x) is the marginal probability density of x, p(ylx) is the conditional 
probability density of y given x, and p(y) is the marginal probability density of y 

P( y)= | P ( y I x ) p( x )dx (A2) 

s 

where the integration is over the entire range S of x, that is, over both the entire 
surface of the earth anti the possible values of frequency. 

The observed range rates are modelled as the sum of the true range rates 
f(x) and a noise component v whose covariance is the matrix Q (see equation IK 
It follows that the conditional probability of y given x is [13] 


p(y| x ) = -== exp . y - f(x )] T O' 1 [y- f ( X )| (A3) 

(2n)“ 2 \ 0 2 



Substituting (A3) into (Al) 


P [ x i n A I y ] = l -— I p(x)cxp -I [ y- f (x) ] T Q' 1 [y- f (x)] dx (A4) 

( 277 )“ 2 y | Q | P(y) \ 2 

The exponential function in (A4) has two maxima. Assume that the area A 
encompasses only one of these, a maximum at x = x A . An approximate expres- 
sion for the integral in (A4) can be obtained by expanding f(x) linearly about 
x A . Let 

x = x A + A x (A5) 

f(x)=f(x A ) + AAx (A6) 

The exponent in (A4) then becomes 

- y [y - f(x A ) - AAx] 1 Q* 1 [y - f(x A ) - AAx] 

= - -y y - f(x A )] T Q _1 [y- f (x A )] - ^ A x T A T Q _1 AAx (A 7) 

The cross product terms 

A x T A T Q _1 [y - f ( x A ) ] + [y - f(x A )j T Q _1 t Ax 
are absent from the right hand side of (A7) because the conditi . that x A is a 
maximum of the exponential function is equivalent to requiring that these terms 
vanish. 

Also, the quantities y - f(x ) appearing in (A 7) can be recognized as the 

A 

residuals It A of the weighted least squares fit to y of f(x A ). Assuming that p(x) 
is a slowly varying function, (A4) becomes 
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1 t -Id 

P [ x i n A I y ] = ^ x A ^ 2 ^a f exp - ^ Ax'A T G -1 AAxdAx (A 8) 

c 2 T7 > M 2 yToT p ( y ) *a 

The integral In (A8) can be approximated by extending the limits of integration 
to infinity, giving for (A 8) 

P(x A ) 1 

P [ x in A I y] = exp-^RjQ _1 R A (A9) 

1 0|^ 2 | A T Q" 1 A | 1 2 p( y ) 

The probability p[x in B/y] is given by a similar expression, with B re- 
placing A throughout. Because of the symmetry of the two maxima and x b , 
the coefficients of the exponential factors are equal, whence equation (36) is ob- 
tained in the form 

exp - ^ RlQ-’R, 

p[x in A I y ] = * (A 10) 

exp - \ r IQ' ,R a + exp - \ RjQ^Re 

by dividing (A‘J) by 1 = P [£ in Aly] + P[x in Bly] . 


Time 

(minutes) 


Table I 

Simulated Data, 10 Hz Resolution 


Noise 

m/s 

Residuals 
of False Min. 
(m/s) 

Residuals 
of True Min. 
(m/s) 

Residual 

Difference 

(m/s) 

0.605 

1.137 

0.704 

0.433 

0.728 

-3.120 

-1.925 

-1.195 

2.859 

3.303 

2.039 

1.264 

-5.109 

-1.321 

-0.818 

-0.503 


Time 

(minutes) 


Table II 

Simulated Data, 1 Hz Resolution 


Range 

Noise 

m/s 

Residuals 

Residuals 

Residual 

Rate 

of False Min. 

of True Min. 

Difference 

m/s 

(m/s) 

(m/s) 

(m/s) 

-1550.382 

-0.225 

0.476 

0.039 

0.437 

0 

-0.189 

-0.303 

-0.105 

-1.198 

1547.107 

0.021 

1.377 

0.111 

1.266 

2876.796 

-0.299 

-0.550 

-0.045 

-0.505 
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